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Abstract. The matrix-valued spherical functions for the pair {K X K,K), K = SU(2), are studied. By 
restriction to the subgroup A the matrix-valued spherical functions are diagonal. For suitable set of represen- 
tations we take these diagonals into a matrix-valued function, which are the full spherical functions. Their 
orthogonality is a consequence of the Schur orthogonality relations. From the full spherical functions we ob- 
tain matrix-valued orthogonal polynomials of arbitrary size, and they satisfy a three-term recurrence relation 
which follows by considering tensor product decompositions. An explicit expression for the weight and the 
complete block-diagonalization of the matrix-valued orthogonal polynomials is obtained. From the explicit ex- 
pression we obtain right-hand sided differential operators of first and second order for which the matrix-valued 
orthogonal polynomials are eigenfunctions. We study the low-dimensional cases explicitly, and for these cases 
additional results, such as the Rodrigues' formula and being eigenfunctions to first order differential-difference 
and second order differential operators, are obtained. 



1. Introduction 

The connection between special functions and representation theory of Lie groups is a very fruitful one, see 
e.g. |3^, For the special case of the group SU(2) we know that the matrix elements of the irreducible 

finite-dimensional representations are explicitly expressible in terms of Jacobi polynomials, and in this way 
many of the properties of the Jacobi polynomials can be obtained from the group theoretic interpretation. In 
particular, the spherical functions with respect to the subgroup S(U(1) x U(l)) are the Legendre polynomials, 
and using this interpretation one obtains product formula, addition formula, integral formula, etc. for the 
Legendre polynomials, see e.g. |2^], 33| for more information on spherical functions. 



In the development of spherical functions for a symmetric pair (G, K) the emphasis has been on spherical 
functions with respect to one-dimensional representations of K, and in particular the trivial representation of 
K. Godement |Q considered the case of higher-dimensional representations of K, see also |3^] for the 
general theory. Examples studied are ||[, @, jl^, 0. However, the focus is usually not on obtaining 
explicit expressions for the matrix- valued spherical functions, see Section ^ for the definition, except for |16 



and In [|16[ the matrix-valued spherical functions are studied for the case {U,K) — (SU(3), U(2)), and 
the calculations revolve around the study of the algebra of differential operators for which these matrix- valued 
orthogonal polynomials are eigenfunctions. The approach in this paper is different. 

In our case the paper p4| by Koornwinder is relevant. Koornwinder studies the case of the compact 
symmetric pair {U, K) = (SU(2) x SU(2), SU(2)) where the subgroup is diagonally embedded, and he calculates 
explicitly vector- valued orthogonal polynomials. The goal of this paper is to study this example in more detail 
and to study the matrix-valued orthogonal polynomials arising from this example. The spherical functions in 
this case are the characters of SU(2), which are the Chebychev polynomials of the second kind corresponding 
to the Weyl character formula. So the matrix- valued orthogonal polynomials can be considered as analogues 



of the Chebychev polynomials. Koornwinder 24 1 introduces the vector- valued orthogonal polynomials which 
coincide with rows in the matrix of the matrix- valued orthogonal polynomials in this paper. We provide some 
of Koornwinder's results with new proofs. The matrix-valued spherical functions can be given explicitly in 
terms of the Clebsch-Gordan coefficients, or 3 — j-symbols, of SU(2). Moreover, we find many more properties 
of these matrix-valued orthogonal polynomials. In particular, we give an explicit expression for the weight, 
i.e. the matrix-valued orthogonality measure, in terms of Chebychev polynomials by using an expansion in 
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terms of spherical functions of the matrix elements and explicit knowledge of Clebsch-Gordan coefficient. This 
gives some strange identities for sums of hypergeometric functions in Appendix Another important result 
is the explicit three-term recurrence relation which is obtained by considering tensor product decompositions. 
Also, using the explicit expression for the weight function we can obtain differential operators for which these 
matrix-valued orthogonal polynomials are eigenfunctions. 

Matrix- valued orthogonal polynomials arose in the work of Krei n [p7| , and have been studied from an 
analytic point of view by Duran and others, see [Q, |^], [|lO|, and references given there. As far as 
we know, the matrix-valued orthogonal polynomials that we obtain have not been considered before. Also 



2 X 2-matrix- valued orthogonal polynomials occur in the approach of the non-commutative oscillator, see |23 
for more references. A group theoretic interpretation of this oscillator in general seems to be lacking. 

The results of this paper can be generalized in various ways. First of all, the approach can be generalized 
to pairs ([/, K) with (Uq)^ abelian, but this is rather restrictive Given a pair ([/, K) and a representation 
6 oi K such that [k\K : ^] < 1 for all representations tt of G and S\m is multiplicity free, we can perform 
the same construction to get matrix-valued orthogonal polynomials. Needless to say, in general it might be 
difficult to be able to give an explicit expression of the weight function. Another option is to generalize to 
{K X K) to obtain matrix-valued orthogonal polynomials generalizing Weyl's character formula for other 
root systems, see e.g. l20| ]. 

We now discuss the contents of the paper. In Section |2] we introduce the matrix- valued spherical functions 
for this pair taking values in the matrices of size (2^+1) x (2^-1-1), £ G ^N. In Section^ we prove the recurrence 
relation for the matrix-valued spherical functions using a tensor product decomposition. This result gives us 
the opportunity to introduce polynomials, and this coincides with results of Koornwinder [ p4[ . In Section 
^ we introduce the full spherical functions on the subgroup A, corresponding to the Cartan decomposition 
U = KAK, by putting the restriction to A of the matrix- valued spherical function into a suitable matrix. 
In Section ^ we discuss the explicit form and the symmetries of the weight. Moreover, we calculate the 
commutant explicitly and this gives rise to a decomposition of the full spherical functions, the matrix-valued 
orthogonal polynomials and the weight function in a 2 x 2-block diagonal matrix, which cannot be reduced 
further. After a brief review of generalities of matrix- valued orthogonal polynomials in Section ^, we discuss 
the even and odd-dimensional cases separately. In the even dimensional case an interesting relation between 
the two blocks occur. In Section ^ we discuss the right hand sided differential operators, and we show that 
the matrix-valued orthogonal polynomials associated to the full spherical function are eigenfunctions to a 
first order differential operator as well as to a second order differential operator. Section ^ discusses explicit 
low-dimensional examples, and gives some additional information such as the Rodrigues' formula for these 
matrix-valued orthogonal polynomials and more differential operators. Finally, in the appendices we give 
somewhat more technical proofs of two results. 
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is supported by Katholieke Universiteit Leuven research grant OT/08/33 and by Belgian Interuniversity At- 
traction Pole P06/02. 



2. Spherical Functions of the pair (SU(2) x SU(2),diag) 
Let K = SU(2), U = K X K and C [/ the diagonal subgroup. An element in K is of the form 

(2.1) fc(a, /3) = ( ^ ) , + = 1, a, /3 G C. 

Let rrit :— fc(e'*/^, 0) and let T d K he the subgroup consisting of the Wj. T is the (standard) maximal torus 
of K. The subgroup T x T C ?7 is a maximal torus of U. Define 

A* = {(mt, m_t) : < < < An} and M = {{mt,mt) : < t < An}. 

We write at = {mt,m-t) and bt = {mt,mt). We have M — Zk,{A^) and the decomposition U = K^^,A^,K^,. 
Note that M is the standard maximal torus of K^, . 
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Figure 1. Plot of the parametrization of the pairs (^1,^2) that contain t upon restriction. 



The equivalence classes of the unitary irreducible representations of K are paramatrized by if — ^N. An 
element £ S determines the space 

:=C[x,y]2£, 

the space of homogeneous polynomials of degree 21 in the variables x and y. We view this space as a subspace 
of the function space C {€? , C) and as such, K acts naturally on it via 

k : p H- !■ p o fc*, 

where fc* is the transposed. Let 

(2.2) j ^ ~£, ~£ + 1, ...,£- 1 J 

We stipulate that this is an orthonormal basis with respect to a Hermitian inner product that is linear in the 
first variable. The representation : K ^ GL{H^) is irreducible and unitary. 

The equivalence classes of the unitary irreducible representations of U are paramatrized by [/ ^ K x K ~ 
iN X iN. An element (^1,^2) G x iN gives rise to the Hilbert space H^^-^^^ := H^^ (K) H'^^ and in turn to 
the irreducible unitary representation on this space, given by the outer tensor product 

T^^'^^fci, fc2)(<^ ^ = r^H^i)(<^) ® T^Hfc2)(<^). 

The restriction of (T^'^'^'^ , H^'^'^^) to K^, decomposes multiplicity free in summands of type £ e |N with 

(2.3) 14 -^2! < ^ < £1 +^2 and £1 + £2 - £ e Z. 

Con versely, the representations of U that contain a given £ E are the pairs {£i,£2) G x that satisfy 
(2.3). We have pictured this parametrization in Figure |l| for £ = 3/2. 



The following theorem is standard, see |25 



Theorem 2.1. The space H^^'^^ has a basis 



: £ satisfies ^ and \j\ < £} 
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such that for every £ that the map (3^^'^^ : — > H^'^'^^ defined by ipj i~> (p^^j^^ is a K -intertwiner. The base 
change with respect to the standard basis {ipjl (X" i^jl} of H^^'^'^ is given by 

where the Cj^'j^'j are the Clebsch-Gordan coefficients, normalized in the standard way. The Clebsch-Gordan 
coefficient satisfies C'^ljlj = if ji + j2 7^ j- 

Definition 2.2 (Spherical Function). Fix a K-type I G and let {ii,£2) G x be a representation 
that contains I upon restriction to K^. The spherical function of type £ G associated to {£i,£2) G x 
is defined by 

(2.4) : U ^ End(FO : ^ ^ {P'/''")* ° T'^'^H^) o 

If is a spherical function of type £ then it satisfies the following properties: 

(1) j, (e) = /, where e is the identity element in the group U and / is the identity transformation in 

(2) <^[^i,^{kixk2) = T'^{ki)<^i{x)T'^{k2) for aU fci,/c2 G and xeU, 

(3) $^ ^ (a;)<I>£ ^ (y) = X£(A;^^)$^ ^ {xky)dk, for all x,y ^ U. Here denotes the character of 
and'x^ = (2^+l)^£. 

Remark 2.1. Definition is not the definition of a spherical function given by Godement Gangolli 
and Varadarajan or Tirao |^ but it follows from property (3) that it is equivalent in this situation. The 
point where our definition differs is essentially that we choose one space, namely End(iJ^), in which all the 
spherical functions take their values, instead of different endomorphism rings for every {/-representation. We 
can do this because of the multiplicity free splitting of the irreducible representations. 

Proposition 2.3. Let EndM{H'') be the the algebra of elements Y G End(ifO such that T^{m)Y = YT^{m) 
for all m E AI . Then g^i^*) C EndM(^f^)- The restriction of to A^, is diagonalizable. 

Proof. This is observation ||2^, (2.6)]. Another proof, similar to Prop. 5.11], uses ma — am for all a G A* 
and m G M so that by (2) 

The second statement follows from the fact that the restriction of any irreducible representation of K^, ^ SU(2) 
to M ~ U(l) decomposes multiplicity free. □ 



The standard weight basis ( |2.2| ) is a weight basis in which ^^^^^U. is diagonal. The restricted spherical 
functions are given by 

(2-5) K.2K)),- E E {c^^^ . 

jl=-il j2=-^2 



which follows from Definition 2.2 and Theorem E.l 



3. Recurrence Relation for the Spherical Functions 

A zonal spherical function is a spherical function for the trivial if-type £ — 0. We have a diffeo- 

morphism U/K^ ~> K : (fci, k2)K^, 1— > /cifc^^ and the left X^-action on U / K.j, corresponds to the action of K 
on itself by conjugation. The zonal spherical functions are the characters on K which are parametrized 
by pairs (£1,^2) with £1 = £2 and we write ipi = Note that ipn = (-l)"^+'(2^ + iy^/'^U2i{cost) by 

1,1,0 



(2.5) and Cj'z'^ ^ = (—1) -'^'(2£+ 1) ^Z^, where Un is the Chebyshev polynomial of the second kind of degree 



n. The zonal spherical function 1^91 plays an important role and we denote it hy ip — (pi. Any other zonal 
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spherical function (p„ can be expressed as a polynomial in cp, see e.g. |p4 , |34[| . For the spherical functions we 
obtain a similar result. Namely, the product of and a spherical function of type i can be written as a linear 
combination of at most four spherical functions of type 

Proposition 3.1. We have as functions on U 



(3.1) 



mi = |£i-i| m2 = |f2-5l 



{■mi,m2),l 



where the coefficients a[^^'^^-* ^ „ are given by 



(mi ,m2)J 



3i,h/ ii,i2,0 ji,ii,ni j2,i2,n2 ni,n2,t 



Jl,j2,'»l,«2,ni,«2 



where the sum is taken over 
(3.3) \ji\<h, 



\j2\<l2, Ki|<i \^2\<^ 



\ni\ < mi and \n2\ < m2. 



2' ' ' - 2 

Moreover, a(^j!^^/2 ^2+1/2) f ^ ^' ^'^^^ ^^^^ ■'"'^"^ (3-^) is a double sum because of Theorem 2.1. 

Proposition 3.1 should be compared to Theorem 5.2 of |2^, where a similar calculation is given for the case 
(SU(3),U(2)). 

Proof. On the one hand the representation T^^'^^ ® T^^h can be written as a sum of at most 4 irreducible 
t/- representations that contain the representation upon restriction to K^. On the other hand we can find a 
'natural' copy of in the space H^^'^^i^H'^'^ that is invariant under the if^-action. Projection onto this 
space transfers via a, defined below, to a linear combination of projections on the spaces in the irreducible 
summands. The coefficients can be calculated in terms of Clebsch-Gordan coefficients and these in turn give 
rise to the recurrence relation. The details are as follows. 

Consider the [/-representation T^i^^2 j^2'2 in the space H^^'^'^ ® H2^2 , By Theorem 2.1 we have 

mi = Ki-i| m2 = K2-^| 

which is a [/-intertwiner given by 

^1 + ^ mi ^2 + ^ 

^€2 ,i ,m2 



- E E E E ^n:C:cl 

mi = |£i-i| ni=-mi ^^^i^^.i | n2=-m2 



Let C H^'^ ® i/2'2 be the space that is spanned by the vectors 



The element ^ 



maps to 



«2 + * 



mi+m2 



ii=-fi j2=- 



■i2=-|mi = |£i-i|»i 



1 1 ni— — m 



1 m2 = |f2-5l "2 = -m2 p=|mi-m2| u=-p 



yM^l ,£2 ,^/^ 2 ' 2 ) 2 '"^1 <^^2 , 2 '^2 ,m2 ,p imi ,1712 

ii,J2,i il,i2,0^ii,ii,ni ^j2,i2,n2 ^ni,n2,u 'Pp,u 



Note that it = ni +712 = ji + ii + j2 + 12 = j, so the last sum can be omitted. Also, since a is a X^-intertwiner, 
we must have p — £. For every pair (mi, 777,2) we have a projection 



j(mi,m2) 



i/™- 

mi=|fi-i| m2 = |£2-^| 



iJ™- 

mi = |£i-i| m2 = K2-5l 
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onto the ^-isotypical summand in the summand H"^^ . Hence 



mi m2 



EV^ Y^ Y^ p£i,£2,^(-i5.5:0 „^l,i,mi „t:i,2'"'^(^n 

2^ 2-^ 2-^ 2-^ 2-^ jlj2,j ■il;i2fi^jl,il,ni -12,^2 ni,"2j '^ij 

ji=-ii 32=^^2 ii=-5 «2=-5 rii=-mi n2=-m2 

The map o a is a ii'^-intertwiner so Schur's lemma implies that 

/I « 1 1 



,m2 ^mi,m2,^ jmi ,m2 



mi m2 



EY^ Y^ Y^ Y^ Y^ ^^i,£2/(-<3.3.0 (-.^i.^.™i^^2,3,m2^mi 

2^ 2-^ 2-^ 2-^ 2-^ jlj2,j il,i2,0^jl,il,ni ^j2,i2,"2 "1, 



m2,£ 

"2 j' 



Jl=-£l j2=-<'2 ii=-5 i2=-5 «l=-mi n2 = -m2 

,(^1.^2) 



is independent of j. Hence it is equal to a 



(mi,m2),^' 



taking j — i. We have 



(^1,^2) ;mi,m2 



^y^^i^j '-^■J>i'o)~ X! "(mi,m2)j'^^J 
mi ,m2 

Moreover, the map 

£1+1/2 £2+1/2 £1+1 £2+^ 

mi = |£i-l/2| m2 = |£2-l/2| mi =|£i - i | m2 =|£2- 5 I 

is a if*-intertwiner. To show that it is not the trivial map we note that aJ^^^^' , l^ is non-zero. Indeed, the 
equalities 



/^*;i,*;2,*; 

Slj2,£ 



iei + £2 + e + iy.A{ei,£2,i) 



(£1- 31)1(^2 -£ + jiy. 



1/2 



,2-^1+2 



^1 + ii + 1 

2£i + 1 



1/2 



and C^''^f+^ 



^1 - Ji + 1 



-,1/2 



2^1 + 1 



where A{£i,£2,£) is a positive function, can be found in Ch. 8] and plugging these into the formula for 
O'li^'l^i \ « shows that it is the sum of positive numbers, hence it is non-zero. 

We conclude that P is non-trivial, so its restriction to is an isomorphism and it intertwines the K^,- 
action. It maps iiT^-isotypical summands to X^-isotypical summands. Hence a\Jff^ = P\Jif^ . Define 



This is a ff-intertwiner. It follows that 

(3-4) ao^f^'^^^ Y: «(™),£/?£"^''"^- 

mi ,m2 

Define the End(if^)-valued function 

^'£i,£2 • U ^ End{H') : x ^ (7^''')* o (r^^-^i ^ t"^-^){x) o 



Note that ^'^^^^^(x) = ip{x)'^[^^i,^{x). On the other hand we have T^^'^" (g) T5^5 = a o 

= E l«£i!l),£p(/5£^''^)*°T^"^'"^°/5£"^''"^- 



ymi,m2 o Q,^ 



Together with (3.4) this yields 

*1 



Hence the result. 



□ 



In Figure g we have depicted the representations (2:2) ^'^'^ (|'3) with black nodes. The tensor product 
decomposes into the four types (| ± i, 3 ± i) which are indicated with the white nodes. 
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Figure 2. Plot of how the tensor product T^^^ (gjT^'i sphts into irreducible summands. 



Corollary 3.2. Given a spherical function there exist 2£+l elements qf[^ g^, j G {~^, —£+!,...,£} ir, 

such that 



(3.5) 



1e[,e2'^{i+j)/2,{i-j)/2- 



The degree of q\'^ is + £2 ~ ^• 

Proof. We prove this by induction on £1+^2- If + ^2 = ^ then the statement is true with the polynomials 

t,j — A 

%-fc)/2,(£+fe)/2 ~ "jy 



-^'■^ = (5j,fe. Suppose £1 + £2 > ^ and that the statement holds for (£'1, 4) with £ < £i + < £1 + ^2- 

inr -J- I (^1-1/2-^2-1/2) 12^ 

We can write |a;/ / ' V^iiM as 



by means of Proposition 3J. The result follows from the induction hypothesis and a|^^ el)^e^^ ^^^^ ^ ^ 



Remark 3.1. The fact that these functions q£[^ £^ are polynomials in cos(t) has also been shown by Koorn- 
winder in Theorem 3.4 of [p4[ using different methods. 



4. Restricted Spherical Functions 
For the restricted spherical functions : A* — > End(i?^) we define a pairing. 



(4.1) 



(^,^)a. = -tr 



$(a) {"^{a))* \D,{a)\da 



where D^{at) — sin^(t), see In Prop. 2.2] it is shown that on the following orthogonality relations 
hold for the restricted spherical functions. 
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Figure 3. Another parametrization of the pairs {ii,£2) containing £; the steps of the ladder 
are paraniatrized by d, the position on a given step by k. 

Proposition 4.1. The spherical functions on U of type i, when restricted to A^., are orthogonal with respect 
to (i-l). In fact, we have 



- + 1)(24 + xY'^''-'-^''- 



(4-2) = IKF^-TTTK^T-^^^^^^XM 



This is a direct consequence of the Schur orthogonahty relations and the integral formula corresponding to 
the U — K^,A^,K^, decomposition. 

The parametrization of the [/-types that contain a fixed K-type £ is given by (2.3). For later purposes we 
reparamatrize (2J) by the function C : N x {-i, ...,£} -J> x iN given by 

'd + e + k d + e-k' 



This new parametrization is pictured in Figure For each degree d have 2£ + 1 spherical functions. By 
Proposition 2.3 the restricted spherical functions take their values in the vector space EndM(-ff^) which is 
2£ + 1-dimensional. The appearance of the spherical functions in 2£ + 1-tuples gives rise to the following 
definition. 

Definition 4.2. Fix a K-type I € and a degree d G N. The function $^ : A* — 7> End(i?^) is defined by 
associating to each point a £ A^ a matrix $^(a) whose j-th row is the vector $^|.^^-|(a). More precise, we have 

(4.3) = ('^'ad,rt(«))^^ ^'^ « € A,. 
The function $^ is called the full spherical function of type i and degree d. 

Let A'^ be the open subset {at ^ A^, : t ^ ttZ}. This is the regular part of A^,. The following Proposition is 
shown in p^ , Prop. 3.2]. In Proposition 5.7 we prove this result independently for general points in A^ in a 
different way. 

Proposition 4.3. The full spherical function $q of type £ and degree has the property that its restriction 
to A'^ is invertible. 

Definition 4.4. Fix a K-type £ € and a degree d € N. Define the function 

(4.4) Q', : a: ^ End{H') : a ^ $^,(a)(<I>^„(a))-i. 
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The j-th row is denoted by Q((dj)(^)- Qd called the full spherical polynomial of type £ and degree d. 
The functions and Q^^^j j.^ are polynomials because (Q^)^ 



is d which justifies the name we have given these functions in Definition 4.4 



q'^^^Mp). The degree of each row of 



We shall show that the functions <f>^ and satisfy orthogonality relations that come from (O). We start 
with the <i>^. This function encodes 2^+1 restricted spherical functions and to capture the orthogonality 
relations of (4.1) we need a matrix valued inner product. 

Definition 4.5. Let he Find{H^) -valued functions on A^,. Define 



(4.5) 



TT 



$(a)(*(a))* \D,{a)\da. 



A, 



Proposition 4.6. The pairing defined by (^.' ) is a matrix-valued inner product. The functions $^ with d G N 
form an orthogonal family with respect to this inner product. 

Proof. The pairing satisfies all the linearity conditions of a matrix valued inner product. Moreover we have 
$(a) (i'(a))* |-D*(a)| > for all a <E A.,. If ($, $) = then $$* = from which it follows that $ = 0. Hence 
the pairing is an inner product. The orthogonality follows from the formula 

(2^+1)2 



{d + e + p+l){d + £-p+l) 



□ 



and Proposition 4.1 
Define 

(4.6) V\a)^¥o{a){¥o{a)yma)l 

with D^(at) = sin^ t. This is a weight matrix and we have the following corollary. 

Corollary 4.7. Let Q and R be Find{H^) -valued functions on A^, and define the matrix valued paring with 
respect to the weight by 



(4.7) 



Q{a)V^{a) {R{a))*da. 



This pairing is a matrix valued inner product and the functions form an orthogonal family for this inner 
product. 



to 



The functions $^ and being defined, we can now transfer the recurrence relations of Proposition |3^ 
these functions. Let Eij be the elementary matrix with zeros everywhere except for the (j,j)-th spot, where 
it has a one. If we write E^ j with \i\ > £ or \j\ > £ then we mean the zero matrix. 

Theorem 4.8. Fix £ G hN and define the matrices A^jB^ and Cd by 



(4.8) 



For a e A* we have 
(4.9) 

and similarly 
(4.10) 



Ad 



Bd = 



Cd = 



E 

e 

E 

e 

E 



C{d,k). 



^C(d-i,fe),^l 



fe,fe+i 



ad.k) 

'■C(d,k-l)J 



E, 



fe,fe-l 



fc=- 



^(a) • $^,(a) = Ad^d+ii^) + Bd'^iic^) + Q$ti(a) 



(p{a) ■ Qiia) = AdQi+i(a) + BdQi(a) + CdQd-M. 
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Note Ad e GL2£+i(K). 

Proo f. It is clear that ( 4.10| ) follows from ( |4.9| ) by multiplying on the rig ht with the inverse of ^q. To prove 
(4.9) we look at the rows. Let p E {—i, — ^ + 1, . . . , ^} and multiply ( [4.9[ ) on the left by Ep^p to pick out the 
p-th row. The left hand side gives (p{a)Ep^p^^^{a) while the right hand side gives 



C(rf,p) 



-,C(d,p) 



Now observe that these are equal by Proposition 3.1 and (4.3). This proves the result since p is arbitrary. □ 



Finally we discuss some symmetries of the full spherical functions. The Cartan involution corresponding 
to the pair {U, K^) is the map 0{ki, = (A:2, fci). The representation T^i-^^ y^2/i q are equivalent via 
the map ® '4'% H> -0^^ ® '^n' follows that 6*<^\^ = This has the following effect on the fuU 



spherical functions from Definition 4.2 
(4.11) 



where J G End(7?^) is given by tp^ i— The Weyl group yV(C/, i^*) = {l,s} consists of the identity and 
the reflection s in G a*. The group 'W{U,K^,) acts on A^, and on the functions on A^, by pull-back. 

J<f>£^ £2 ("^^ •^'^'^ '^^^ a e A* . The effect on the full spherical functions of 



Lemma 4.9. We have s*<i>^^ ^^(a) 
type £ is 

(4.12) s*^'d = '^dJ- 

Proof This follows from (fj) and the fact that C-^j^j = (_i)'?i+«2 

Proposition 4.10. The functions <i>^ commute with J. 



-jl,-j2,-j' 



□ 



Proof. The action of 6* and Sa on A* is just taking the inverse. Formulas (4.11) and (4.12) now yield the 
result. □ 



5. The Weight Matrix 



We study the weight function : A^ End{H^) defined in (16), in particular its symmetries and explicit 
expressions for its matrix elements. First note that is real valued. Indeed, commutes with J, 



(5.1) JV'ia)J=j¥oia)Jij¥oia)jr\D4a)\=¥oia)¥oiar\D,{a)\ = V'ia), 
since J* = J and = 1. This also shows that V is real valued, 

(5.2) V^^=V'ia-t)^JV'iat)J = V'iat). 



Lemma 5.1. The weight has the symmetries Vp ^ = V^p = V^p 



Vlg^^p for p,qe {-£,...,£}. 



Proof. The first equallity follows since V (a) is self adjoint by (4^) and real valued by (5^). Since V (a) 



commutes with J we see that Vp^^ — V_p _g. 
Set 

(5.3) v'ia,) = ^i{a)¥,{ar 

so that V^{at) = v^{at)\D^,{at)\ — v^{at) sin^ t. Note that for ~t <p,q < t the matrix coefficient 

(5.4) v'-{at)p.q = tr {•^\_+p {at) f ^'1+, {at)) 

\ 2 ' 2 V 2 ' 2 / 

is a linear combination of zonal spherical functions by the following lemma. 



□ 
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Figure 4. Plot of the decomposition of the tensor product T^^^ ^ T'^-i into irreducible 
representations. The big nodes indicate the irreducible summands, the big black nodes the 
ones that contain the trivial iiT^-type upon restricting. 



Lemma 5.2. The function U — )■ C : a; i— ?► tr £2^^^ (*^mi,m2 (^)) *) bi-K -invariant function and 

min(^i+mi ,£2+^12) 

(5-5) tr($^^_^Jat)($i,,™,(at))*) = Yl c„C/2„(cost) 

n— max(|^i — mi | , 1^2 — ^2 | ) 

if ii + TOi - {£2 + e Z and tr {^^[^^i,2^at) ('I'Li,m2 ("«))*) = otherwise. 

Proof. It follows from Property (2) that the function is bi-_R'- invariant, so it is natural to expand the function 
in terms of the zonal spherical functions U2n corresponding to the spherical representations T"'", n € ^N. 
Since T^^ '^^ is equivalent to its contragradient representation, we see that the only spherical functions occurring 
in the expansion of tr {^\^ 12^^) m2(^))*) ^^"^ ones for which {n,n) ^ A — {(71,1,71,2) £ |N x : 
£i+mi~ni G Z, \£i — mi \ < ni < ^1 + mi, ^ l£ ^2 < £2 + ^2} since the right hand side corresponds to 
the tensorproduct decomposition T^i'^^ (g) r™!^'"^ = 0(ni,n2)eA ^"1^"=^ , see Proposition |^ and Figure |. □ 

Given d, e G N and —£ <p,q<£ wc write (^{d,p) = (^1,^2), Ci^Tl) = (wijW2)- Then we have 

where the sum is taken over 

|j|<^, \ji\<h, \J2\<£2, Nil < mi, h\<m2, 
satisfying ji + ^2 = ii + ?2 = j- This equals 

(5.6) (<i>^,(a,)«(a,))*)^^= ^ 4^,^1.-2,36''^* 

|s|<min(£i+mi,£2+m2) 

where 

/r 7\ _ \ " f ^ii,e2/^rni,m2,A'^ 
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where the sum is taken over 



\j\<i, \jl\<il, |J2|<^2, Nil < mi, K2|<TO2, 

satisiying ji + j2 — ii + 12 ~ j and j2 — ji + ii — 12 — s. Since [/„(cost) = e~™* + e~*("^^)* 
follows from (5.6) and Lemma |3.2| that we have the following summation result. 

*«i,£2,mi,m2,s independent of s. 



it 



Corollary 5.3. Let \s\ < ma.x{\£i — mi|, {£2 — ?Ti2|). T/ien 



Note that the sum in (5/7) is a double sum of four Clebsch-Gordan coefficients, which in general are 3F2- 
series 

We now turn to the case £1 = ^-^^,£2 = fni — ^^,^2 = Because of Lemma 5J. the next theorem 
gives an explicit expression for the weight matrix. 

Theorem 5.4. Let q — p < and q + p < 0. For n — 0, ...,£ + q there are coefficients cfi{p, q) G Q>o such 
that 

(5.8) = sin2(i) J2 cUp, q)U2i+p+q-2n{cos{t)) 

where the coefficients are given by 

2£ + l {£-q)l{£ + q)\ {p~£)i+g. 



(5.9) 



£ + P + l 



{£ + p + 2)e+<j-, 



By Cn+p+q{p, q) ~ Cn{~q, —p) the expansion ( pT^ ) with (5.9) remains valid for q < p. 

Proof. Since min(£i + mi, £2 + ^2) — £ + and max(|^i — mi\, \£2 — 'Ti2|) = in this case we find the 
expansion of the form as stated in (5.8). It remains to calculate the coefficients. Specializing (5.6) and writing 



(5.10) 

we find 
(5.11) 



(cs; 



+m)/2,{i-m)/2,£ 



^3,31 



+32 ' 



( i+m \ ( £~m \ 
\]i + (e+m)/2) \j2 + {e~m)/2) 



(A) 



E 



VpqicOSt) 

with 

(5.12) 4(p,q) 



i+Ep- / min(^±a,r+^±£) 

J2 4(p,9)exp(z(-2(z+j»)= { E ' 



Fl^_^{p,q) e 



-2irt 



r=-(£+£±a) V=max(- 



£ + p 



j + {£ + p)/2j yi + i£ + q)/2 



£ + q 



min(-j+-2E,j+^) 

E 

k—msix{—j— .i~^-2^ ) 



k-j+{e-p)/2 



) {k-l+U-q)/2) 



( Y 



From this we can obtain the explicit expression of v^{at)p^q in Chebyshev polynomials. The details are 
presented in Appendix A. □ 

Proposition 5.5. The commutant 

{V\a) : a £ A.,}' {Y e End(ij'^) : V'^{a)Y ^ YV^{a) for all a e A^} ^ {v^{a) : a £ A,}', 
is spanned by the matrices L and J . 



Proof. By Proposition 4.10 we have J S {V^{a) : a G A*}'. It suffices to show that the commutant contains 
no other elements than those spanned by / and J. 



Let v^{at) = X]n=o ^"■('^°^^)^"' ''^ith An € Mat2£+i(C) by Theorem bA Then for B in the commut ant 
it is necessary and sufficient that AnB — BAn for all n. First, put C — A2e- Then by Theorem 5.4 



Cp^g = (^^^p) ^(5p,_g. The equation BC = CB leads to Bp^g = {C-^BC)p^ 



1 ~ C„.-„^-P--9 



C Q , —qC. 



9. — q'-^-Qig 



Pi — P*— ' — P iP 
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by iteration. Since Cq.-q ~ Cp-p if and only ii p — q ov p — —q we find i?p,g — for p ^ q ov p ^ —q. 



Moreover, Bp,p — B-p,^p and Bp,-p — B-p,p. 



Secondly, put C = ^2£-i- Then, by Theorem 5.4, we have 



SO the non-zero entries are different up to the symmetries C'^ ^ = CLp^^g = C'^p = C'_g _p. Now BC = C'B 
implies by the previous result Bp ^Cp ^ + Bp^-pC'_p q = C'p^Bq^q + C'p _^Bq^q. Take q — 1 — p to find Bp^p ~ 
Bi^p^i^p = -Bj,-i,p-i unless p = or p = f , and take q = p — 1 to find Bp .^p — i?i-p,p-i — unless 
p = or p = f . In particular, for ^ e 5 + this proves the result. In case £ G N we obtain one more equation: 
Bq o — Bi i + Si.-i. This shows that B is in the span of / and J. □ 

The matrix J has eigenvalues ±1 and two eigenspaces and H^. The dimensions are [£ + 1/2 J and 
\£ + 1/2] . A choice of (ordered) bases of the eigenspaces is given by 

(5.13) {^Pj -i^ij ■■ -e<j <0,£- j eZ} and {V^j + : < j < ^, £ - j € Z}. 



Let Yg be the matrix whose columns are the normalized basis vectors of ( ^.13 ). Conjugating with Yg yields 



a matrix with two blocks, one block of size \£ +1/2] x \£ + 1/2] and one of size [£ + 1/2] x [£ + 1/2] . 

Corollary 5.6. The family (Q^)d>o (ind the weight are conjugate to a family and a weight in block form. 
More precisely 

Y, Q,{at)Y, - Q ),Y, V (^at)Y, - ^ ^^^^^^ 

The families (Q^ -|-)d>o orthogonal with respect to the weight Vj.. Moreover, there is no further reduction 
possible. 

Proof. The functions can be conjugated by Yg. Since the commute with J we see that the Yf^Q^^i 
has the same block structure as Yf'^V^Yg. The blocks of the Yf'^Q'^Yi are orthogonal with respect to the 
corresponding block of Y^'^V^Yi. The polynomials _ take their values in the (— l)-eigenspace Hi_ of J, the 
polynomials in the (+l)-eigenspace of J . The dimensions are + 1/2] and \£ + 1/2] respectively. 
A further reduction would require a n el ement in the commutant {V^{a) : a € A*}' not in the span of / and 



J. This is not possible by Proposition 5.5. □ 



The entries of the weight with the Chebyshev polynomials of the highest degree 2£ occur only on the 



antidiagonal by Theorem 5.4. This shows that the determinant of v (at) is a polynomial in cost of degree 
2£(2£+ 1) with leading coefficient (—1)^(2^+1) Hp^-^ co(p, —p)2^^ ^ 0. Hence is invertible on away from 
the zeros of its determinant, of which there are only finitely many. We have proved the following proposition 



which should be compared to Proposition 4.3 



Proposition 5.7. The full spherical function <I>g is invertible on A^, except for a finite set. 



In particular, is well-defined in Definition 4.2, except for a finite set. Since is polynomial, it is 
well-defined on A. 

Mathematica calculations lead to the following conjecture. 

Conjecture 5.8. det(i;^(at)) = (1 - cos^ i)^(2^+i) np=_^(22^c^(p, -p)). 



Conjecture ^^is supported by Koornwinder [gj. Prop. 3.2], see Proposition O, Conjecture 5^ has been 
verified for £ < 16. 
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6. The matrix orthogonal polynomials associated to (SU(2) x SU(2), diag) 
The main goal of Sections |3[ ^ and |^ was to study the properties of the matrix valued spherical functions of 



any K-type associated to the pair (SU(2) x SU(2), diag). These functions, introduced in Definition 2.2, are the 



building blocks of the full spherical functions described in Definition 4.2. We have exploited the fact that the 
spherical functions diagonalize when restricted to the subgroup A. This allows us to identify each spherical 
function with a row vector and arrange them in a square matrix. 

The goal of this section is to translate the properties of the full spherical functions obtained in the previous 
sections at the group level to the corresponding family of matrix valued orthogonal polynomials. 

6.1. Matrix valued orthogonal polynomials. Let be a complex x matrix valued integrable function 
on the interval (a, b) such that W is positive definite almost everywhere and with finite moments of all orders. 
Let Matjv(C) be the algebra of all iV x iV complex matrices. The algebra over C of all polynomials in the 
indeterminate x with coefficients in MatAr(C) denoted by Mat^v (C) [x] . Let (•, •) be the following Hermitian 
sesquilinear form in the linear space Mat jv (C) [x] : 

(6.1) {P,Q) ^ I P{x)W{x)Q{x)*dx. 



The following properties are satisfied: 

(1) {aP + bQ,R) = a{P,R) +b{Q,R), for all P,Q,Re MatA,(C)[x], a, 6 e C, 

(2) (TP, Q) = T{P, Q), for aU P, Q e MatAr(C)[x], T e MatAr(C), 

(3) (P,Q)* (Q,P), for all P,Q e Mat^v (C) [a;] , 

(4) (P, P) > for all P e Mat at (C) [a;] . Moreover if (P, P) = then P = 0. 

Given a weight matrix W one constructs a sequence of matrix valued orthogonal polynomials, that is a 
sequence {Rn}n>o, where P„ is a polynomial of degree n with nonsingular leading coefficient and (P„, Rm) = 
if n ^ m. 

It is worth noting that there exists a unique sequence of monic orthogonal polynomials {Pn}n>o in 
MatAr(C)[a;]. Any other sequence of {Pri}ri>o of orthogonal polynomials in MatAr(C)[a;] is of the form 
Rnix) = AnPn{x) for some An e GLAr(C). 

By following a well-known argument, see for instance one shows that the monic orthogonal 

polynomials {P„}„>o satisfy a three-term recurrence relation 

xPn{x) = P„+i(x) -f Bn{x)Pn{x) + C„P„_i(a;), n > 0, 

where Q-i = and P„, C„ are matrices depending on n and not on x. 

There is a notion of similarity between two weight matrices that was pointed out in Q. The weights W 
and W are said to be similar if there exists a nonsingular matrix M, which does not depend on x, such that 
W{x) = MW{x)M* for aU x e (a, b). 

Proposition 6.1. Let {Pri,i}n>o be a sequence of orthogonal polynomials with respect to W and M £ GLn{C). 
Then the sequence {Rn,2{x) = Rn.i{x)M^^}n>Q is orthogonal with respect to W = MWAI*. Moreover, 
if {P,i.i} *s the sequence of monic orthogonal polynomials orthogonal with respect to W then {Pn,2{x) = 
MPn^i{x)M^^} is the sequence of monic orthogonal polynomials with respect to W. 

Proof. It follows directly by observing that 

Rn,2{x)W{x)R„^^2{x)*dx = [ P„4 (a;)M- ^ W?(a;) (M" ^ ) * P,„,i (x) * dx 



Rn.i{x)W {x)Rm,i{x)* dx — 0, if n ^ m. 

The second statement follows by looking at the leading coefficient of P„^2 and the unicity of the sequence of 
monic orthogonal polynomials with respect to W. □ 



MATRIX VALUED ORTHOGONAL POLYNOMIALS 



15 



A weight matrix W reduces to a smaller size if there exists a matrix M such that 

Wix) ^ M (^'l^'^^ ^r^^^)^M*, foraUxe (a,6), 

where Wi and W2 are matrix weights of smaller size. In this case the monic polynomials {Pn}n>a with respect 
to the weight W are given by 



Pnix) = M 



n D I \ ] ' " > 0, 



where {P„,i}„>o and {Pn,2\n>a are the monic orthogonal polynomials with respect to Wi and W2 respectively. 

6.2. Polynomials associated to SU(2) x SU(2). In the rest of the paper we will be concerned with the 
properties of the matrix orthogonal polynomials Qd- For this purpose we find convenient to introduce a new 
labeling in the rows and columns of the weig ht V. More precisely for any £ e let W be the (2£ + 1) x (2£ + 1) 
matrix given by 



(6.2) ^Jl- x'^ W{x)n,m = V{a^^ccosx)-l+n-l+r,i, n, 771 G {0, 1, . . . , 2£}. 

It then follows from Theorem ^.4| that 

if n < m and W{x)n,m = W{x)m^n otherwise. 

We also consider the sequence of monic polynomials {Pd}d>o given by 

(6.4) -Pd(a;)„,„i = T^^Q<j(aarccosa;)-f+n -f+m, n, TO G {0, 1 , . . . , 2f } , 



where is the leading coefficient of the polynomial Qd{0'a,iccosx)i which is non-singular by Theorem 4.8. Now 
we can rewrite the results on Section |^ in terms of the weight W and the polynomials Pd- 

Corollary 6.2. The sequence of matrix polynomials {Pd{xy\d>o is orthogonal with respect to the matrix valued 
inner product 

{P,Q)^[ P{x)W{x)Q{x)*dx. 



Theorem 4.8 states that there is a three term recurrence relation defining the matrix polynomials Qd- These 
polynomials are functions on the group A. We can use ( |6.4[ ) to derive a three term recurrence relation for the 
polynomials Pd. 

Corollary 6.3. For any £ G the matrix valued orthogonal polynomials Pd, are defined by the following 
three term recurrence relation 

(6.5) xPd{x) = Pd+i{x) + T^'BdTdPdix) + T^'CdTd-iPd-iix), 

where the matrices Ad, Bd and Cd are given in Theorem J^.S and taking into account the relabeling as in the 
beginning of this subsection. 

6.3. Symmetries of the weight and the matrix polynomials. In this section we shall use the symmetries 
satisfied by the full spherical functions to derive symmetry properties for the matrix weight W and the 
polynomials Pd. 

For any n G N, let /„ be the n x n identity matrix and let J„ and Fn be the following n x n matrices 

TL—l n— 1 

(6.6) Jn^Y. = J2{-iyE,,,. 

i=0 4=0 

For any nx n matrix X the transpose X* is defined by {X*)ij = Xji (reflection in the diagonal) and we define 
the reflection in the antidiagonal by {X'^)ij — Xn-j.n-i- Note that taking transpose and taking antidiagonal 
transpose commute, and that 

(XY = {X-^Y = X"^' = JnXJn. 
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Moreover, {XZY = Z^X''- for arbitrary matrices X and Z. We also need to consider the {2£ + 1) x {2£ + 1) 
matrix Y defined by 



(6.7) 




if ^ = 



2n+l 



if 



Proposition 6.4. The weight matrix W{x) satisfies the following symmetries 

(1) W{xY = W{x) and W{x)'^ = W{x) for all x e [-1, 1]. Thus 

J2l+lW{x)J2i+l = W{x), 

for all X Cz [—1,1]. 

(2) W{-x) = F2i+iW{x)F2e+i for all x e [-1,1]. 
Here F21+1 is the {21 + 1) x [21 + 1) matrix given in (6.6). 



Proof. The symmetry properties oiW in (1) follow directly from Lemma 5.1 



The proof of (2) follows from (6.3) by using the fact that Un{~-x) = (— l)"[/„(a;) for any Chebyshev 
polynomial of the second kind Un{x), so that W{—x)n,m — (— l)"'''"'VF„^m(a;). □ 



The weight matrix W can be conjugated into a 2 x 2 block diagonal matrix. In Corollary 5.6 we have 
pointed out this phenomenon for the weight V. The following theorem translates Corollary 5^ to the weight 
matrix W. 



Theorem 6.5. For any I £ ^N, the matrix W satisfies 

W{x) = YW{x)Y^ = 



WAx) 






W2{x) 



where Y is matrix given by (3.7). Moreover if {Pd,i}d>Q (resp. {Pd.2}d>o) is a sequence of monic matrix 
orthogonal polynomials with respect to the weight Wi{x) (resp. W2{x)), then 

'Pd.i{x) 

Pd.2, 



Pd{x) 



d>0, 



is a sequence of matrix orthogonal polynomials with respect to W . There is no further reduction. 

The case £ — {2n + l)/2, n G N, le ads to weights matrices W of even size. In this case W splits into two 
blocks of size £ + In Corollary 6.6 we prove that these two blocks are equivalent, hence the corresponding 



matrix orthogonal polynomials are equivalent. 

It follows from Proposition 6.4 (1) that there exist (n + 1) x {n 
A{xY = A{x) and 

' A{x) B{x) 



1) matrices A{x) and B{x) such that 



(6.9) 

for all X e [-1, 1]. 
Corollary 6.6. Let 



W{x) = 

(2n+ l)/2, n e Z. Then 
YW{x)Y* = 



dt 



B{xY' A{x) 



Wi{x) 
W2{x) 



where 

Wi{x) = A{x) + B(a;) J„+i, W2{x) = J„+iF„+iH/i(-a;)F„+i J„+i. 
Here A{x) and B{x) are the matrices described in ( |6.3[ ) and (6^). Moreover, if {Pd,i}d>o is the sequence of 
monic orthogonal polynomials with respect to Wi{x) then 

(6.10) Pd,2{x) = {~iyJn+lFn+lPdA~x)F„+iJ„+i, 
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is the sequence of monic orthogonal polynomials with respect to W2 (x) . 

Proof. In this proof we will drop the subindex in the matrices Jn+i, Fn+i and we will use J and F instead. 
It is a straightforward calculation that for {n + I) x (n + l)-matrices A, B C and D the following holds 

(A B\ t _ I fA + + J(C + B"**) B-C'^' + {D'^' - A)j\ 
\C d) ^ 2 \J{D'^*' ~ A) + C - B'^* D + A'^* ~ + B'^*)j) 

In particular for the weight function W we get 



YWf-rWt -Y ( ^^^'^ B{x)\^t _ (A(x)+B{x)J 
ivv(x)i - y^^^yt A{x)'iy ^\ J{A{x) ~ B{x)J)J 

This proves that 

Wi{x) = A{x) + B{x)J, W2{x) = J{A{x) - B{x)J)J. 
It follows from Proposition |6.4| (2) that A{—x) — FA{x)F and B{—x) — FB{x)F. Therefore we have 

JFWi{-x)FJ = JFA{-x)FJ + FJB{-x)FJ = JA{x)J - JB{x) = W2{x). 
This proves the first assertion of the theorem. 



The last statement follows from Proposition 6.1 □ 



7. Matrix valued differential operators 

In the study of matrix- valued orthogonal polynomials an important ingredient is the study of differential 
operators which have these matrix- valued orthogonal polynomials as eigenfunctions. In this section we discuss 
some of the differential operators that have the matrix- valued orthogonal polynomials of the previous section 



as eigenfunctions. The calculations rest on the explicit form of the weight function (6.3). 



7.1. Symmetric differential operators. We consider right hand side differential operators 
(7.1) D^Y.d^F,{x), d^-^, 

in such a way that the action of D on the polynomial P{x) is 

s 

PD = Y,d\P){x)F^{x). 

i=0 

In |l^, Propositions 2.6 and 2.7] one can find a proof of the following proposition. 

Proposition 7.1. Let W — W{x) be a weight matrix of size N and let {Pn},i>o be the sequence of monic 



orthogonal polynomials in MatAr(C)[a;]. If D is a right hand side ordinary differential operator as in (7.1) of 
order s such that 

PnD = AnPn, foralln>0, 

with An G A, then 

i 

F, = F,{x) = ^^h ^3 e Matjv(C), 
is a polynomial of degree less than or equal to i. Moreover D is determined by the sequence {A„}„>o and 



An = ^N.^'(^), foralln>0, 



i=0 



where [n\i = n{n — 1) • • • (71 — i + 1), [n]o = 1. 
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We consider the following algebra of right hand side differential operators with coefficients in MatAf(C)[x]. 

V^{D ^^d'F, : F, e MatA,(C)[x], degi^, < i). 

i 

Given any sequence of matrix valued orthogonal polynomials {Rn}n>n with respect to W , we define 

V{W) = {DeV : RnD = Tn{D)Rn, Tn{D) G MatAr(C), for all n > 0}. 
We observe that the definition of T>{W) does not depend on the sequence of orthogonal polynomials {Rn}n>o- 

Remark 7.1. The mapping D i-> r„(D) is a representation of T>{W) in for each n > 0. Moreover the 
family of representations {r„}„>o separates the points oi'D{W). Note that I?(M^) is an algebra. 

Definition 7.2. A differential operator D d V is said to be symmetric if {PD, Q) ~ {P, QD) for all P,QCz 
MatAr(C)[x]. 

Proposition 7.3 i^). If D eV is symmetric then D € V(W). 

The main theorem in says that for any D £1) there exists a unique differential operator D* e 'D{W), 
the adjoint of D, such that {PD,Q) = {P,QD*) for ah P,Q e MatA,(C)[x]. The map D ^ D* is a 
operation in the algebra 'D{W) . Moreover we have 'D{W) = S{W) ® iS{W), where S{W) denotes the set of 
all symmetric operators. Therefore it suffices, in order to determine all the algebra 'D{W), to determine the 
symmetric operators S{W). 

The condition of symmetry in Definition 7.2 can be translated into a set of differential equations involving 
the weight W and the coefficients of the differential operator D. For differential operators of order 2 this was 
proven in Theorem 3.1]. 

Theorem 7.4. Let W{x) be a weight matrix supported on (a, 6). Let D £ V be the differential operator 

D = d^F2ix)+dFi{x)+F°, 
Then D is symmetric with respect to W if and only if 

(7.2) F2W = WF2, 

(7.3) 2{F2Wy = WF^ + FiW, 

(7.4) {F2W)" - {FiW)' + Ff)W ^ WFq, 
with the boundary conditions 

(7.5) lim F2{x)W{x) = 0, lim {F2{x)W{x)y - Fi{x)W{x) = 0. 

7.2. Matrix valued differential operators for the polynomials P„. As in the previous section, we will 
denote by {P„}fi>o the sequence of monic orthogonal polynomials with respect to the weight matrix W . We 
can write the weight as W{x) — p{x)Z{x) where p{x) = (1 — a;)^ (f + a;)^ and Z{x) is the (2£ + 1) x (2£ + 1) 
matrix whose (n, m)-entry is given by 

_ (2£+l)(2£-m)!m! " (n-2^(2£ + 2-^ 

(7.6) Z(a:)„,„ - — -j-^ — ^(-1) ^ 2)™_, t\ 

= c(n, m, t)Un+m~2t{x), 

t=Q 

if n < m and Z{x)n,m = Z{x)m,n otherwise. 

Once we have an explicit expression for the weight matrix W we can use the symmetry equations in Theorem 
7.4 to find symmetric differential operators. If we start with a generic second order differential operator 

2 i 

D^Y. P^i^) = ^'^h e Matw(C), 

j=0 j=0 
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then the equations (7^), (7^) and {TA^) lead to hnear equations in the coefficients Fj. It is easy to solve these 
equations for small values of N using any software tool such as Maple. We have used the general expressions for 
small values of N to make an ansatz for the expressions of a first order and a second order differential operator. 
Then we prove that these operators are symmetric for all N by showing that they satisfy the conditions in 
Theorem 



7.4 



In the following theorem we show the matrix polynomials Pn satisfy a matrix valued first order differential 
equation. This, phenomenon which does not appear in the scalar, case has been recently studied in the 
literature (see for instance [||, Q). 

Theorem 7.5. Let E be the first order matrix valued differential operator 

d 



E = 

where the matrices Ai [x) and are given by 

21 



dx 



Ai{x)+Ao, 



i=0 
2i 



2e^ 



21 



E. 



i=0 



21 



Ei,, 



-^i.i—l 1 



Then E is symmetric with respect to the weight W ; hence E € D{W). Moreover for every integer n>0, 

Pn{x)E ^ An{E)Pn{x), 

where 

n{e-i) {2e + 2){i -2ey 



21 



21 



Ei, 



Proof. The proof of the theorem is performed by showing that the differential operator E is symmetric with 
respect to the weight W . It follows from Theorem 



7.4 



with F2 — 0, that E is symmetric if and only if 



(7.7) 
(7.8) 

with the boundary condition 
(7.9) 



W{x)Ai{x)* + Ai{x)W{x) = 0, 
-(Ai(x)VK(x))' + A^W{x) = W{x)Al, 



lim Ai{x)W{x) 



The second statement will then follow from Propositions 7.1 and 7.3, 

The verification of (7.7) and (7.8) involves elaborate computations, see Appendix ^ 

Theorem 7.6. Let D be the second order matrix valued differential operator 



□ 



D = {l-x')^ + 



dx^ 



Bi{x) + Bo, 



where the matrices Bi{x) and Bq are given by 



21 



B,{x)=Y. 



i=Q 
21 



(4^ + 3)(i-2^) 
21 



^ {2£ + 3){i-2e) 



y-^ [i - 2ti{it - 2P - U - i) 

Bq — — E!i: 



i=0 



21 



Then D is symmetric with respect to the weight W{x); hence D € D{W). Moreover for every integer n>Q, 

Pn{x)D = A„{D)Pn{x), 
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^n{D) = 2^ n(n - 1) + — Er, 



whe 

Proof. The proof of the theorem is similar to that of Theorem 7.5, see Appendix ^ □ 
Corollary 7.7. The differential operators D and E commute. 

Proof. To sec that D and E commute it is enough to verify that the corresponding eigenvalues commute. The 
eigenvalues commute because they are diagonal matrices. □ 



As we pointed out in Theorem 3.5, for any £ E the matrix weight W and the polynomials P„ are 
{2£ + 1) X {2£ + 1) matrices that can be conjugated into 2x2 block matrices. More precisely 

YW(x)Y - I ^ W2{x) h ^^"^^^^ -I P,M 



where Y is the orthogonal matrix introduced in ( |6.7D and Wi , Wi are the square matrices described in Corollary 
6.6. Here {Pn,i}n>o and {Pn,2}n>o are the sequences of monic orthogonal polynomials with respect to the 
weights W\ and respectively. 

Proposition 7.8. Suppose £ = {2n + l)/2 for some integer n, then E splits in (n + 1) x [n + 1) blocks in the 
following way 

where 

El = (^^^Ii(x)+Io, 

E2 ~ Fn+lJn+lAi{ — x)Jn+lFn+l + J„+i^o J„+i F„+i . 

Here Fn+i, Jn+i o,fs the matrices introduced in ( |6.6| ). The matrices Ai and Aq are given by 

'^i2£-i)^ V-(^ V- « z:. (2<'+1)j:, 
Ai{x) - - 2^ — E,^n-'i.-l + X ^ J Ei^n-i + —E,^n^,+i H — E^^n+l, 

1=0 i=0 4=1 

^0 — 7, J^i,n-i- 

i=0 ^ 

Proof. The proposition follows by a straightforward computation. □ 
Proposition 7.9. Suppose £ £ N, then we have 



yFeAi{x)Fi V2 OmJ \FiAo{x)Fe vq -{£+l)hj 




where Ai and Aq are n x n matrices given by 

i=0 1=0 i=l 



j=0 
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and the vectors vo,vi,V2 & are 
wo = (0,0,--- ,0), vi 



,0,---,0 , V2=\- 



{2£+l)V2 
4(^+1) ' 



O,--- ,0 



Proof. The proposition follows by a straightforward computation. 



□ 



Let us assume that £ = {2n + l)/2 for some n G N so that the weight W and the polynomials P„ are 
matrices of even dimension. Proposition says that 



(7.10) 

A simple computation shows that 



P„{x)E = YAnY'Pn{x), n>0. 



yA„F* = -(£+!)/: 



2)1+2 



A„,i 
A„,2 



where A„.i is a (n + 1) x (n + 1) matrix (depending on n) and 

An, 2 — Fn+lJn+l-A7i,lJn+lFn+l- 



It follows from (7.10) that the following matrix equation is satisfied 



-{£+l)Pn^l{x) PnAx)Ei 
PnAx)E2 -{£+l)PnAx) 



-{£+l)Pn,l{^) K.lPnAx) 



Therefore the polynomials Pn,i and Pn,2 satisfy the following differential equations 



(7.11) 
(7.12) 



Pn,lEi — A„ iP„ 2 — 0, 
Pn,2E2 — A„,2^ri4 = 0. 



Finally it follows from ( 7.11 ), ( 7.12 ) and ( 6.10 ) that for every n > 0, the polynomial is a solution of the 
following second-order matrix valued differential equation 

PnsEiE2 — A„,iA„,2^ri,l — 0. 

We can also obtain a second order differential equation for Pn,2- 



8. Examples 

The purpose of this section is to study the properties of the monic orthogonal polynomials {Pn}n>o presented 
in Section |^ for small dimension. For £ = 0, ^, 1, |, 2, we show that these polynomials are solutions of certain 
matrix valued differential equations. We will show that the polynomials can be defined by means of Rodrigues' 
formulas and we will give explicit expressions for the three term recurrence relations. 



8.1. The case £ = 0; the scalar weight. In this case the polynomials {Pn}n>o are scalar valued. The 
weight W reduces to the real function 

Wix) ^ {l~x)^l + x)^, xe[-l,l]. 

Therefore the polynomials Pn are a multiple of the Chebyshev polynomials of the second kind: P„(x) = 
2-"C/„(.t), n e N. 
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8.2. The case £ — weight of dimension 2. In this case the polynomials {Pn}ri>o are 2x2 matrices. 
The weight W is given by 



It is a straightforward computation that 



3 



^ ' V {l~x)2{l + x)2j V2 V-1 1 

Observe that Wi{x) — {I — x)i {I + x)i and W2{x) — (1 — a;)^ (1 + a;)^ are Jacobi weights and therefore we 
have 

2"n!(?i + 2)! (13),, „ 2"n!(n + 2)! (a i) 

where {Pi"'''''}n>o are the classical Jacobi polynomials 



8.2.1. Differential equations. By Theorem 7.5 we have 

AP,.W(- l,)+P„w(-o=' °)^(-V' ,1)^ 

We can conjugate the differential operator E by the matrix Y to obtain 
The monic polynomials 



dx \x-l / 2 V 1 -1/ ' 



satisfy 

-3/2 n + 3/2 



Pn{x)E = A„f„(x), where A„(x) = . ^ _3/2 

Now the fact that is an eigenfunction of E is equivalent to the following relations between Jacobi 

polynomials 

(1 + + - (n + ^)Pi^'^\a;) = 0, 

da; 2 2 

(1 - x)^Pii^'^'\x) + ^Pi^'^^(x) - (n + l)pi^-'i\x) = 0. 
da; 2 2 

8.3. Case t — 1; weight of dimension 3. Here we consider the simplest example of nontrivial matrix 
orthogonal polynomials for the weight W . The weight matrix W of size 3 x 3 is obtained by setting 1=1. 
We have 

J 1 / 3 3^ 4a;2 - 1^ 

(8.1) VF(a;) = (l-a:)2(l + x)2 3a; + 2 3a; 

\4a;2-l 3a; 3 



We know from Theorem |6.5| that the weight W{x) splits into a block of size 2x2 and a block of size 1x1, 
namely 

/w/ i \ n \ + 2 3V2x 

YW{x)Y' = r^^}^''^ , J = (l-a;)5(l + x)H 3^/2:E a;^ + 2 

V V 4(1 -.2) 

From Theorem |6.5| the monic orthogonal polynomials Pn{x) with respect to W{x) reduce to 

'PnA{x) 



P-n. 



P„,2(a;) 
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where {Pn,2}n>o are the monic polynomials with respect to Wi{x) and {Pn.2}n>o are the monic polynomials 
with respect to the weight W2{x). 

Remark 8.1. The weight W2 is a multiple of the Jacobi weight (1 — x)"{l + x)^ corresponding to a = 3/2 
and /3 = 3/2. The monic polynomials {P„,2}n>o are then a multiple of the Gegenbauer polynomials 

_ 2"n!(n + 3)! (1,1) 
^"'^^''>- (2n + 3)! " ^''>- 

8.3.1. The first order differential operator. By Theorem 7.5 we have that the monic polynomials P„ are 
eigenfunctions of the differential operator E. More precisely the following equation holds 

^ , -a; 1 0\ /-4 0\ - 4 0^ 

— P„(a:) -i i +P„(x) -2 = -2 \ Pn{x). 



dx 



2 2 

-1 x/ \ 0/ \ On, 



Now we can conjugate the differential operator E by the matrix Y to obtain a differential operator i? = F £' F*. 
The fact that the polynomials P„ are eigenfunctions of E says that the polynomials P„ are eigenfunctions of 
i?. In other words 

d~ /O ^\ ^ /-2 2\ /-2 0n + 2\ 

^ +Pn{x) -2 = -2 P„(x). 
"^"^ \x -V2 J \2 0-2/ V" + 2 -2/ 

We can now rewrite the equation above in terms of the polynomials P„_i and P„,2- 

^Pn,2{^){x -a/2) +P„,2(a;)(2 0) = (n + 2 0)P„,i(x). 

Since Pn,2 is a Gegenbauer polynomial, we see that the elements of the first row of P„,2 can be written explicitly 
in terms of Gegenbauer polynomials. 

8.3.2. Second order differential operators. In this subsection we describe a set of linearly independent differ- 
ential operators that have the polynomials P„.i as eigenfunctions. 

Proposition 8.1. The matrix orthogonal polynomials {Pn, i}„>o satisfy 

PnsD, = A„p,)P„,i, J = 1, 2, 3, n > 0, 
where the differential operators Dj are 

n , d"^ \ fx^ -2x\ f d \ fhx -&\ /4 

T'2 



dx^J \ % -I J \dxj \ l J \0 

_(f_\ f -2x 8x^ ~ 4\ f d_\ f-8 S2x\ fO 16 
dx^ [x'-2 2x r\dx [Gx -4j"*'U 



and the eigenvalues A, are given by 



A ^ >('^ + 4) + l \ /(n + 2)2 



3) 



8(n + 2)(n + l)' 

^{n + 3){n + 2) 

Moreover, the differential operators Di, D2 and satisfy 

DiD2^D2Di, DiDs^DsDi, D2D3 ^ D3D2. 
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Proof. The proposition follows by proving that the differential operators Dj, j — 1,2,3 are symmetric with 
respect to the weight mat rix W i. This is accomplished by a straightf orw ard computation, showing that the 
differential equations (7.2), (7.3), (7.4) and the boundary conditions (7.5) are satisfied. As a consequence of 
Remark 7.1, the commutativity properties of the differential operators follow by observing the commutativity 
of the corresponding eigenvalues. □ 

8.3.3. Rodrigues' Formula. 

Proposition 8.2. The matrix orthogonal polynomials {Pn,i{^)}n>o satisfy the Rodrigues' formula 



i.2) 



where 



Pn,l{x) = C 



(l-x2)2+" 



4a;2 + 2 'SV2x 
3\/2x + 2 



2n 
n+2 
_ ^/2nx 

71+1 



x/2n 



n+2 
_ n 
ri+1/ 



in) 



_ (-l)"2-2«-2(„ + 2)(n + 3)0F 
^~ (2n + 3)r(n + f) ' 

Proof. The proposition can be proven in a similar way to Theorem 3.1 of [lOl. We include a sketch of the 



proof for the for the sake of completeness. First of all, we recall that the classical Jacobi polynomials P„ 
satisfy the Rodrigues' formula 

('—1 
2"n! 



>{x) 



Let R{x) and y„ be the matrix polynomials of degree 2 and 1 



R{x) 



+ 2 3V2x 
3V2x + 2] ' 



Ynix) 



2n 

n+2 n+2 
_ \/2nx _ n 



n+1 



n+\ , 



so that the (S^) can be rewritten as 



Pn,\{x) = C 



(l-a;2)2+» (i?(x) + r„(.T)) 



(») 



i^r'(^)- 



Then by applying the Leibniz rule on the right hand side of (5.2), it is not difficult to prove that 

P„,i(x) = c\n(n - 1)[(1 - + a;)5+"]("-2)(l _ + x)-'^ R!' [x)R{x)-^ 

+ cn[(l - a;)5+"(l + a;)H»]("-i) (i _ a;)-5(l + x)-^(ii"(x) + r'(a;))ii'(x)-i 

+ c[(l - x)5+"(l + a;)^+"](")(l - x)-^(l + a;)-^(/ + Y{x)R{x)-^). 

Now by applying the Rodrigues' formula for the Jacobi polynomials we obtain 

(8.3) P„,i(a;) = 2"n!(-l)"c[pi^'^^ (a;)(/ + y (a;)i?(x)-i) 

- i(l - x)(l + x)P<y '(x)(y'(a;) + i?'(x))i?(a:)-i 

+ i(l-x)2(l + x)2p(if (x)i?"(x)i?(x)-i]. 

Now with a careful computation we can show that the expression above is a matrix polynomial of degree n 
with nonsingular leading term. Using integration by parts it is easy to show the orthogonality of Pn,i and x™, 
m = 0, 1, . . . , n — 1, with respect to the weight W\. □ 
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8.3.4. Three term recurrence relations. In Corollary 3.3 we show that the matrix polynomials Pn{x) of any 
size satisfy a three term recurrence relation. The recurrence relation for the polynomials Pn.i can then be 
obtained by conjugating the recurrence relation for Pn(x) by the matrix Y. The recurrence coefficients 
are given in terms of Clebsch-Gordan coefficients and are difficult to manipulate. For ^ = I we can use the 
Rodrigues' formula (3.2) to derive explicit formulas for the three term recurrence relation for the polynomials 

Pn,l- 

First we need to compute the norm of Pn,i{x). The Rodrigues' formula (3.2) and integration by parts lead 



to 



p II 2 o^2ri— 1 




8(n+l)^ 



If {Pn,i\n>o is a sequence of orthonormal polynomials with respect to Wi with leading coefiicients r2„, then 
it follows directly from the orthogonality relations for the monic polynomials that = ri~^(il*)~^. The 

orthonormal polynomials 7^„,i with leading coefficient 








v^(ra+3) 



satisfy the three term recurrence relation 
where A„ = fln^iilz:^ and 



The coefficient of : 



Bn = f2„[coef. of x" ^ in P„ i ~ coef. of x" in P„+i i]r2„ 
in Pn^i can by obtained from (^.3|). Now a careful computation shows that 



1 / n{n+3) 

2 V (n+l)(n+2) 









n(n+3) 



Bn — 





2(n+l)(n+4) 



(n+2)(n+3) 



2(n+l)(n+2)^ 

Therefore the monic polynomials P„^i satisfy the three term recurrence relation 

xP„,l{x) Pn+is{x) + BnPn.lix) + CnPn-l.l{x), 

where 







B„ 



2(n+l)(n+2) 



(n+2)(n+3) 





' n(K+3) 
4(n+2)(ri+l) 





4(n+2)^(n+l)'^ 



8.4. Case i = 3/2; weight of dimension 4. The weight matrix W of size 4 x 4 is obtained by setting 
I = 3/2. 

/ 4 4a; |(4x2 - 1) 4x(2a;2 - 1)\ 

4a; |(4a;2 + 5) |x(2a;2 + 7) |(4a;2 - 1) 

4x 
4 



VF(a;) = (l-a;2)2 



|(4a;2-l) |a;(2a;2 + 7) |(4x2 + 5) 
V4a;(2a;2 - 1) |(4a;2 - 1) 



„2 1^ 4^,^2 

We know from Corollary |6.5| that the weight W{x^ splits in two blocks of size 2x2, namely 

'VFi(x) 



W{x) = YW{x)Y^ 







W2{X) 



where 



and 



W,{x) = 4(1 - xfl^X + xf/' i") ' 1 



i(4a;-l) 
(2a;2 + 2a; + 5) / ' 



W2{x) = J2F2Wi{-x)F2J2, where F2 



1 

-1 



and J2 



1 

1 
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It follows from Corollary 3.6 that the monic orthogonal polynomials Pn,2 with respect to the weight W2 are com- 
pletely determined by the the monic orthogonal polynomials Pn,i with respect to Pn,2{x) — J2^2^n,i(— a;)^2 ^2- 
Therefore we only need to study the polynomials P„^i. 

8.4.1. Differential operators. In this subsection we describe a set of linearly independent differential operators 
that have the polynomials Pn.i as eigenfunctions. 

Proposition 8.3. The matrix orthogonal polynomials {P„,i}n>o satisfy 

PnsDj=K^{D,)Pn,u J = l,2,3,n> 0, 
where the differential operators Dj are 

dx^ \dx [-1 6x-2r lo ' 



Di = - 1) 



_cP_\ f -3x + 3 9x^ -9x 
dx^J \x^ + x-2 3a; - 3 

d \ / -9 36a; - 18\ / 18 



D. = 



,dx J \8x + 4 -3 ; V12 

and the eigenvalues Aj are given by 

/ 9(n + 2)(n + ir 

''^^> ~ \in + A){n + 3) 

Moreover, the differential operators Di, D2 and satisfy 

DiD2=D2Di, DiDs^DsDi, D2D3 ^ D^D^. 

8.4.2. Rodrigues' Formula. The monic orthogonal polynomials {Pn,i(a^)}n>o satisfy the Rodrigues' formula 

13 1 

P„,i(a;)=c (l-a;)2+"(l + a;)2+"(i?(x)+r„(a;)) W^\x), 



where 



and 



R{x) 
%,{x) 



2-^"-^(^l)"(n + 3)(n + 4)0F 

2a;2-2a; + l i(4a;-l) 
i(4a;- 1) i(2a;2 + 2a; + 5)y' ' 

;iT3 3(n+3) (22; + 1) 

3(ri+l) (1 ~ 2x) ~3(ji+l) . 



8.4.3. Three term recurrence relations. The orthonormal polynomials Vn,i{x) — ||^ri,i|| ^Pn,i, with leading 
coefRcient 

/ 22n + l(„+l) 

7T(n+4) " 



9(n+l)^22"+i(»+2) 



— 

1 ^ , 

^7r(n+3)(n+4) 

satisfy the three term recurrence relation 



xl^n (,X 
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where 



1 / »(n+4) 

2\J («+l)(n+3) ^ 



7t(w+4) 



2(n+2)^(n+l)(n+3); 







3 



Bn = 



2->/(n+l) (n+2) (n+3) (n+4) 
3 2 



^2^(Ti+l)(n+2)(Ti+3)(n+4) (n+2)(Ti+3) 

Therefore the monic polynomials Pn,i{x) satisfy the three term recurrence relation 

xPn,l = -Pn+1,1 + BnPn,l + CnPn-1,1, 

where 

U Tvr-rWTT^ \ _ / 4(n+l)(n+3) ^ 



\ 2{n+l)(n+2) (n+2)(n+3) / V ^ 



rt^(n+4)'^ 



4(ri+l)(ri+2)2(n+3), 



8.5. Case i = 2; weight of dimension 5. In this subsection we consider the 2x2 irreducible block in the 
case £ = 2, where the matrix weight W is of dimension 5. This case completes the list of all irreducible 2x2 

blocks obtained by conjugating the weight W by the matrix Y. 
The 2x2 block is given by 

Wi{x) = {l-x)Hl + x)^(^ lOx 16x^+4)- 
As before, we denote by {Pn,i}n the sequence of monic orthogonal polynomials with respect to Wi. 

8.5.1. Differential operators. In this subsection we describe a set of linearly independent differential operators 
that have the polynomials P„ j as cigenfunctions. 

Proposition 8.4. The matrix orthogonal polynomials {Pn,i}n>o satisfy 

P^^^Dj = A„(£),)P„,i, j = 1,2,3, n > 0, 
where the differential operators Dj are 







o) 




''5 






S) 








+ ( 



5 



and the eigenvalues Aj are given by 



in{n + 6)-3 \ /(n + l)(n + 5) 



/ ^(n + 5)(n + 4) 

^^(n + 2)(n+l) 

Moreover, the differential operators Di, D2 and D3 satisfy 

D1D2 = D2D1, D1D3 D^Di, D^Ds ^ i?3-D2. 
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8.5.2. Rodrigues' Formula. The monic orthogonal polynomials {Pn,i{x)}n>o satisfy the Rodrigues' formula 



(1 - a;)2+"(l + a;)2+"(i?(x) + r„(a;)) 



in) 



where 



and 



(-l)"2-2"-4(n + 3)(n + 4)(n + 5)0F 
(2n + 5)r(n+ I) ' 



R{x) = 



lOx 



lOx I6x^ 



3n 



7 J^n[X) — I 6na; 12k 
ri+4 ri+4 



8.5.3. Three term recurrence relations. The orthonormal polynomials Vn,iix) — ||i^n,i|| ^Pn,i with leading 
coefRcient 



-)2n+2 \/("+2)("+l) 







^7r(n+4)(ri+5) 

Q g„ V2(k+1) 

•^7r(rt+5)y 



satisfy the three term recurrence relation 

xPji (x 

) = + BnVn{x) + A*7'„_i(a;), 

where 

^(" + 5) Q 



2^(ri+l)(n+2)(n+3)(n+4) 



^(«+l)(n+4)^ 
2 



^(ri+l)(ri+2)(n+4)(n+5) 



^ V("+l)("+2)("+4)(«+5) 

Therefore the monic polynomials P,i,i(a;) satisfy the three term recurrence relation 

xPn,l = Pn+IA + BnPn^l + CnPn-1,1, 

where 





(n+4)(T!+5) 



2(n+l)(Ti+2) n — I 4(n+l)(n+2)(n+3)(n+4) 

j' ^--[ 

Appendix A. Transformation formulas 





^(?^+5) 



4(n+l)(n+4), 



The goal of this appendix it to prove Theorem 5.4. We use the standard notation for the Pochhammer 
symbols and the hypergeometric series from In the manipulations we only need the Chu-Vandermonde 
summation formula ^ Cor. 2.2.3] which reads 



-n, a 



;1 



(c- a)r 



{C)r. 



(A.l) 2F1 

and Sheppard's transformation formula for 3-F2's [y, Cor. 3.3.4] written as 

{-n)k{a)n{h)r 

e -I- K)n-k\<^ -I- l^jn-k 

k=0 

{-n)k{a)k{a + b - n - d - e + l)k 

u — a)n-k\e — Ujn~k 

k=0 



(A.2) 



y^(e + k)n-k(d + k)n-k 

k=0 

- a)n-k{e - a)n-k 
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Proposition A.l. Let £ e and p,q ^ such that \p\,\q\ < i, £ — p,£ — q E 'Z, q — p < and q + p < 0. 
Let s e {0, . . . ^£ + q} and define 



(A.3) 

Then we have 

(A.4) ei{p,q) 



n / \n + s 



E 



2 ■ 



i£ + p+l) {2£y. 



EM) 



T=0 



{i + p + 2)^+,_tT! 



Proof. First we reverse the inner summation using M — £ ^ p — s — m to get 
(A.5) 



£+q-s 
n=0 



+ p\ £ + q 

n J \ n + s 



( 21 \2 
M=0 \i_p_M+n) 



We rewrite the inner summation: 

M J \i-p-M- 



(A.6) 



21 



M=Q \I-p-M+n) 



^ig^(-ir(-f+P). 



2£ 

£ — p + n 



-1 l—p—s 

E 

M=0 



{-£ + p + s)m{£ + P - n + 1)m 



Ml{-£ + p-n) 



M 



B{M) 



where B{M) = {I — p — M + l)n{p — q + M + s + l)e+q-s-n is a polynomial in M of degree i + q — s that 
depends on i,p,q,n and s. The polynomial B{M) has an expansion in (—!)*(— M)t, 



(A.7) 



B{M) = YAf{-in-M), 



The coefficients At = At{£,p, q,n, s) can be found by repeated application of the difference operator Aa// 
f{M + f ) - /(Af). Let Alf be its i-th power. We have 



(A8) 

In other words, 
(A.9) 



B = At 



M=0 



i+q 



Af=0 



We calculate ( [A.6| ) by substituting (A.7) in it. Interchanging summations, the inner sum can be evaluated 
using ([A.lj) (after shifting the summation parameter). We get 



2£ 

£ — p + n 

{~2£ - l)e-p-s-ti-i + p + s)t{e + p-n+l)t 
{-£ + p - n)i_p_s 



i-p^s (i-p\( i-q \ ^, 



/ 21 

M=0 \e-p-M+nl 



E^= 

t=0 
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Substituting ( |A.10| ) and (A^) in (A^) and simplifying gives 
(A.11) 



t\ 



t=Q 

e+q-s 



{p-q + s + M + l)i+g^s J2 



i-e~q- sU-e - pUe -p-m + i\ 

nii-e - p ^ t)„{-£ - p - M)„ 



The inner sum over n is a 3i^2-series, which can be transformed using ( |A.2[ ). Note that the t-order difference 
operator can now be evaluated yielding only one non-zero term in the sum over n. This gives 



. {-£ + p)s+t{.2 + £-q + s + t)e+g-s-t 



£ + p+l {2£y. ^' ^£ + p + 2),+t{£ + q-s-ty. 

Reversing the order of summation using T~l + q — s — t yields 

(2^ + 1) {£ ~ q)\{£ + g)! ,+g_r jp - ^),+,-t(2 + 2£~ T)t 

^ ^ ' {£ + p + 2),+g^TT\ 



(A.13) 



{£ + p+l) {2£)\ 



T=0 



as was to be shown. 



□ 



Proof of Theorem 5.4- We already argued that there is an expansion in Chebyshev polynomials (q^)- From 
(5.11) it follows that there are coefficients df.{p, q) such that 



(A.14) 

The coefficients d^{p,q) and c'^{p,q) are related by 
(A.f5) 



vlg{cost)= d'^{p,q)e- 

r=-(£+£±a) 



diip,q)= ciip,q). 



n=0 



Let q — p< 0, q + p < and r > and substitute r{s) — s + in ( |5.f f] ). Comparing this to (A.f4) shows 



(A.16) 



n=0 ^ / \ ' / „j^(; 



Vm+s/ 



+ S / ^ ( 21 \ 
r?i-0 \m+n+s) 



ioT s — Q,. . . ,£ + q. Now we use Proposition A.l to show that q) equals 
(A.17) 



i+q-{r+3^) 



{2£+l) {£-q)\{£ + qy/- ^ ^ ' ,^^_„ (p - f (2 + 2^-n)„ 

^ ^ ' {£ + p + 2)t+g_^n\ 



{£+p+l) {2£)\ 



n=0 



for r = H^, . . . , ^ + It follows that 



(A.I8) 

This proves the theorem 



U \^ {£~q)\{£ + q)\ {p^£)i+q-n , ,+,_ j2£ + 2-n), 

''"^^''^^ £ + p+l (2£)! {£ + p + 2),+,.r} ' n\ 



□ 



We can reformulate Proposition A.l in terms of hypergeometric series. 
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Corollary A. 2. For N e N, a,b,c e E so that 0<a<N,0<b<N, and additionally a<b, N<a + b 
and 0<c<N — bwe have 



E 



{-C)m {b+ l)m {b+ l)m 

(N - a~ c + 1),„ ml (6 - N)r 



-b,N-a-b-c,N-b-m + l,N-b-m + l 
N — b — c + 1, —b — m, —b — m 



fN+l\ 



E 



{-a)N-b-n (-1)^-"-" (A^ + 2 - n)„ 
{N - a + 2)N-b-nnl 



(w-a-c) ^N-b-c) 71=0 

The 4F3-series in the summand is not balanced. Note that the case s = leads to single sums, and the 4^3 



boils down to a terminating 2F1 which can be summed by the Chu-Vandermonde sum, so Corollary A. 2 
be viewed as an extension of Chu-Vandermonde sum (A.l). 



The coefficients d^{p,q) of (5.11) with \r\ < 2—2 are independent of r. Corollary 5.3 in case i — £1 + I 



mi + m2 can be stated as follows. 

Corollary A.3. For N e N,. a, 6, c G N so thatO < a < N , < b < N , b < a, a + b < N, andO<c< N-a-b 
we have 



N y 

a+cJ 



n=0 



n\ (c + 1)„ [a + c- N)^{a + c- iV)„ 



-a, — a — c, A — a — c — 71 + 1, A — a — c — 71 + 1 
X 4-r3 I ; 1 

A — a — — c+1, —a — c — 71, —a — c — 77 



A + 1 /A 



-1 b 

E 

m— 



(-a)™(A-6 + 777 + 2)f,_„(-ir 



TV-a+lV&y (A-a + 2)„(&-777)! 

In particular, the left hand side is independent of c in the range stated. 



Different proofs of Corollaries A. 2, A. 3 using transformation and summation formulas for hypergeometric 
series have been communicated to us by Mizan Rahman. 



Appendix B. Proof of the symmetry for differential operators 



Proof of Theorem VLq. In terms of p{x) and Z(x), the equations (7.7) and (7.8) are given by 



(B.l) 
(B.2) 



Z{x)Ai{xy + Ai{x)Z{x) 
p{x') 



= -A^{x)Z{x) 



p{x) 



Ai{x)Z{x) - Ai[x)Z'{x) + A^Z{x) - Z{x)Aq. 



As a consequence of the properties of symmetry of the weight W, it suffices to verify the conditions above for 
all the (ti, 77T,)-entries with 77 < ?7i. Here we assume that n < m. The case n = m can be done similarly. The 
first equation (B.l) holds true if and only if 



^n,m— ('^)m.m—l ^~ -^n.rn^l ('^)m.7n ~t" ^n,m+l^l (■^)m,m+l 

+ Zn-l.mAi{x)n.n-l + Zn^rnAi{x)n.n + Zn+l.mAi{x)n.n+l = 0, 
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for all n < m. In order to prove the expression above we replace the coefficients of Ai and Z in the left hand 
side and we obtain 



m— 1 ^ m 

> ^ c(n,m~ l,t)Un+m-2t-i{x) J — 2_^c(n,m,t)xUn+m-2t(xj 



21 



t=0 



t=0 



21 -m 
2i 



771 + 1 



c(n, m + 1, i)C/„+m_2t+i(x) " ^ X! ~ ™' t)Un+7n-2t-l{x 



t=0 



e-n 

£ 



c(n, m, a; C/„+,„-2t( 



t=o t=o 
By using the recurrence relation xUr{x) — ■^Ur^i{x) + ^Ur+i{x) we obtain 



Zt — n 

a;) H 2^ c(n + 1, m, t)Un+m~2t+i\^ 



m— 1 

E 

«=0 



— c(n,m- - ^c(n- l,m,t) 



2i? — n — m 



2^ 



■ c(n, rn, t) 



{x) 



E 

t=0 



2£-TO , ^ , 2^-n , . 2£-n-m , 
c(n, TO + 1, t ) H cyn + 1, to, t ) cln, m, t ) 



n+m-2t+l 



(x) 



n 2£ — m 2£ — m — n 
-^c(n, m + 1, to) — c{n, TO + 1, m + 1) c(n, m, m) 



A simple computation shows that the coefficient of t/n-m-i in the expression above is zero. Now by changing 
the index of summation t we obtain 



(B.3) 



2£-fn , ^ 2£-n , ^ 2£ - n - m , 

— — — c{n, TO + 1, 0) H c(n + 1, to, 0) c(n, m, 0) 



Un+m+l{x) 



E 

t=0 



m n 2^ — n — TO, 

— c(n,m- l,t) - ^c(n- l,TO,t) c{n,m,t) 



2£~m , ^ 2£-n , ^ 2£~n~m , 

— — — c(n, TO + 1, t + 1) H c(n + 1, m, t + 1) to, t + 1) 



n+m-2t-l 



Using the explicit expression of c(n, m, in ( [7.(]| ) we obtain that (B.3) is given by 

(TO+l-f + n)(2^-TO + l) (-n + l + 2£)(TO + l-t + n) 2£~n-m 



t=o 



2£{~m +l + t-n + 2£) 2£{-m + 1 + t - n + 2£) 



+ 



(2£+ 1 -t)(TO+ 1) (2^ + 1 -t)(n + 1) (2£-n-TO)(TO + 1 -f + n)(2£+ 1 -t) 



2£{t + l) 



2£{t + l) 



£(-TO + l + f-?i + 2£)(i+l) 



2£ 



Un+m-2t-l{x) = 0, 



since the sum of the terms in the square brackets is zero. This completes the proof of (B.l). 
Now we prove (B.2). The (n,TO)-entry of the right hand side of (B.2) is given by 



x{Ai{x)Z{x))n,m - (1 - x'^){Ai{x)Z'{x))n^m. 

Using ( |7.6D we obtain 

(B.4) (1 - X^)[{Ao)n.n - A[{x)n.n - {Ao)rn,m]Zn 



+ (1 - x'^)[{Ao)n.n - A[{x)n.,n - (^o)m,m]^n 



E" £(n-TO+l)-TO 2\TT I \ 
c(n,TO,t)(l - a; )(7„+m-2t(2;)- 

t=0 * 
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(B.5) x{Ai{x)Z{x))n^rn. = ^ [~ '^'^^'^ ~ l,m,t)xUn 

e-n 



+m-2t-l 



(x) 



t=0 



2i — n 

^ C(n, m, t) x'^ Un+m~2t{x) H 2£~^^^ ^ ^ Un+m-2t+l{x) 



(B.6) {l-x^){A,{x)Z'{x)) 
£-n 



7n 



2£ — n 

c{n,m,t)x{l - x'^) K+rn-2tix) H — + l,m,t) (1 - x^) U^_^_ 



m-2t+l\ 



Now we proceed as in the proof of the condition (BJ). In (B^) and ( |B.5| ) we use the thre e te rm recurrence 
relation for the Chebychev's polynomials to get rid of the factors x and x^. Equation (B.6) involves the 
derivative of the polynomials U. For this we use the following identity 

Finally we change the index of summation t and we use the explicit expression of the coefficients c(n, m, t) to 
complete the proof. 

□ 



The boundary condition (7.9) can be easily checked 



Proof of Theorem \7.q . We will show that the conditions of symmetry in Theorem 7.4 hold true. The first 
equation (7.2) is satisfied because A2{x) is a scalar matrix. Equation (7.3) can be written in terms of p{x) 
and Z(x) in the following way 

(6x - Bi{x))Z{x) + 2{x'^ - l)Z'{x) - Z{x)Bi{xy = 0. 



This can be checked by a similar computation to that of the proof of Theorem 7.5 



Now we give the proof of the third condition for symmetry. If we take the derivative of (7.3), we multiply 
it by 2 and we add it to ( |7.4| ) we obtain the following equivalent condition 

(B.7) {W{x)Bi{x)* - Bi{x)W{x))' - 2{W{x)Bo - BoW{x)) = 0. 

We shall prove instead that 



W{x)Bi{x)* -Bi{x)W{x) 



2 (^j W{x)dx^ Bo - 2Bo (^j W{x)d3^ = 0, 



which is obtained by integrating (B.7) with respect to x. Then (B.7) will follow by taking the derivative with 
respect to x. 



We assume n < m. The other cases can be proved similarly. We proceed as in the proof of (B.l) in Theorem 
7.5 to show that 



(B.8) {W{x)Biixy - Bi{x)W{x))n^r. 



-p[X) C{n, m,t) nl ^ , ^ , „ , o/;\ I i\ Un+m-2t-l{x) 



£{-m + l + t-n + 2£) {t + I) 



On the other hand we have 

(B.9) - f 2 ( / W{x)dx] Bo + 2Bo[ I W{x)dx 



E 

t=0 



2c{n,m,t){{Bo)r 



{Bo)n,n) I P{x)Un 



%+m-2t 



(x)dx 
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It is easy to show that the foUowing formula for the Chebyshev's polynomials holds 



Therefore we have that ( P3.9D is given by 

m1n^ I i ^^ (■Bo)m,m - (■Bo)«,«) ( c{n, ■m,t+l) \ 

(B.IO) p(x)^c(n,m,0 ^ ^ _ e(n,m,t) " V ^"^-^-i^") 

Now ( [B.IO ) is exactly the negative of ( |B.8| ). This completes the proof of the theorem. □ 
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